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CN \ It is shown that the study of correlations in the associative production of B c and D 

Qh mesons at LHC allows one to obtain essential information about the B c production 

< 

CL, 1 I. INTRODUCTION 

Recent measurements of B c meson mass and lifetime in CDF jl] and DO [2] experiments 

in ■ 

£> ■ are the first steps in the experimental research of quarkonia with open flavor. The measure 
. ment results are in a good agreement with the theoretical predictions for the B c mass 

m CDF = 6>2 756 ± o.0029(stat.) ± 0.0025(sys.) GeV, 
m°° = 6.3000 ± 0.0014(stat.) ± 0.0005(sys.) GeV, 



3H5|: 



m *«* = 6.25 ± 0.03 GeV; 



as well as for the decay time 



rgDF = .448i°;°i(stat.) ± 0.032(sys.) ps, 
r°° = 0.475+°;^(stat.) ± 0.018(sys.) ps, 

^heor = q 4g ± q Q5 pg 



Unfortunately, the experimental estimations of the cross section value were not published. 
Thus, the mechanism of B c meson production can not be understood from the obtained data 
due to poor experimental statistics, as well as due to large uncertainties in the theoretical 
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predictions. Only the planed experimental research at LHC, where about 10 10 events with 
B c mesons per year are expected, will improve the situation. This huge amount of events will 
allow to obtain the information on the production cross section distributions, on the decay 
branching fractions, and in some cases, on the distributions of decay products. Obviously, 
the new experimental data on B c meson from LHC experiments have to stimulate new 
theoretical studies. 

At first sight the production of B c meson can be represented as the heavy b quark produc- 
tion followed by the fragmentation of b quark into B c meson by analogy with the production 
of B mesons, at least for transverse momenta larger than the B c meson mass 1 . But, as it 
was shown in the framework of pQCD, there is a large nonfragmentation contribution into 
this process js-ll]. The nonfragmentation terms in the production amplitude enlarge the 
total cross section and increase the ratio between yields of B* and B c states by a factor of 
two even at large transverse momenta. 

However, our previous analysis showed that it is quite difficult to obtain a clear experi- 
mental evidence of the nonfragmentation contribution into the production |8j. In this article 
we try to fill the gap in our understanding of B c production and show that the study of 
correlations in the associative production of B c and D meson at LHC could be an essential 
information source of B c production mechanism. The cross section distributions of the B c 
meson production in hadronic and gluonic interactions have been obtained by us within the 
numerical calculation method developed in our previous studies js]. 



II. FRAGMENTATION AND RECOMBINATION CONTRIBUTIONS INTO B c 

PRODUCTION 

The B c production amplitude within the discussed approach can be subdivided into 
two parts: the hard production of two heavy quark pairs calculated in the framework of 
perturbative QCD and the soft nonperturbative binding of b and c quarks into quarkonium 
described by nonrelativistic wave function. The production amplitude within this approach 

1 Within fragmentation approach it is assumed that bb phase space can be factorized from the total phase 
space. This assumption is not valid for final state B c + b + c at small momenta, because c quark is massive. 
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can written as follows: 



A SJj * = J T^( Pi , k{q)) ■ ■ C%-±^, (1) 

where T,r s - is an amplitude of the hard production of two heavy quark pairs; 
\1/^ 2 is the quarkonium wave function; 

J and j z are the total angular momentum and its projection on z-axis in the B c rest frame; 

L and l z are the orbital angular momentum of B c meson and its projection on z axis; 

S and s z are B c spin and its projection; 

CA Z are Clebsh-Gordon coefficients; 

Pi are four momenta of B c meson, b quark and c quark; 

q is three momentum of b quark in the B c rest frame (in this frame (0, q) = k(q)). 
Under the assumption of small dependence of T^ s t on k(q) 



(2) 



and particularly for the 5* wave states 

A~R s (0).T blc5 ( Pi )\^ , (3) 

where Rs{0) is a value of radial wave function at origin. 

The calculations within the discussed technique are the most simple for the process of 
B c production in the e + e~ annihilation, where only four leading diagrams contribute to the 
process amplitude (see Fig. [T]). As it was shown in Q, [12I 13|, the diagrams (3) and (4) are 
suppressed as m 2 c /m1 for the S wave state production due to the large gluonic virtuality. 
The special choice of gluonic field gauge allows to neglect the contribution of diagram (2) 
at large energies yfs ^> Mg c . The remained contribution of diagram (1) can be naturally 
interpreted as the b quark production followed by the fragmentation of b quark into B c 
meson. Thus, in e + e~ annihilation at large energies the consideration of leading diagrams 
for the B c meson production leads the well known factorized formulas for the cross section 
distribution over z = 2Eb c / ' \fs and over B c transverse momentum pt'- 

^ = m-D(z), (4) 
1 

da Bc f da bi fp T \ D(z) ^ 



dpT J dkx 
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Figure 1: The leading order diagrams for the 
process e + e~ — >• i? c + 6 + c. 



Figure 2: The fragmentation functions for the 
processes b — ^ B* (solid curve) and b ^ B c 
(dashed curve) calculated within pQCD. 



where is a transverse momentum of 6-quark. 

The analytical forms of fragmentation functions for S wave states are known from 
(see Fig. E]): 
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A-B e (*) 



:(6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )^ 2 



llnml (l-(l-r)zY 

- 2(1 - r)(6 - 19r + 18r 2 )^ 3 + 3(1 - r) 2 (l - 2r + 2r 2 )^ 4 ), (6) 



D- b ->B*Az) 



2a 2 | J R s (0)| 2 rz(l - z) 2 , , 9n „ 

sl (2-2(3-2r)^ + 3(3-2r + 4r 2 )^ 2 



277rm3 (1 - (1 - rb) 



- 2(1 - r)(4 - r + 2r 2 )^ 3 + (1 - r) 2 (3 - 2r + 2r 2 )z 4 ), (7) 



where a s is a strong coupling constant and r = m T ^ c m ■ 

The graphics presented in figures [21121 have been obtained for the following parameter 
values: a s = 0.2, m c = 1.5 GeV, m b = 4.8 GeV, R s (0) = 2.14 GeV 3 . As it is shown in [14 1 
the uncertainties in these parameters leads to the quite large uncertainties in the absolute 
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Figure 3: The leading order diagrams for the process gg — >• B c + b + c. 



cross section values. In this work we study the production characteristics which weakly 
depend on these parameters: the cross section distribution shapes and the ratio between B* 
and B c yields. 

The relative yield of B* and B c in e + e~ annihilation obtained within pQCD calculation 
R^+ e - = cr(B*)/a(B c ) ~ 1.4. Thus the naive spin counting which fairly predicts this ratio 
for B* and B [R + _ ~ 3) can not be applied to B c and B c production. 2 

At first sight it would be reasonable to assume that for the gluonic B c production the 
fragmentation mechanism is also dominant at least from transverse momenta larger than 



B r mass. But as it was shown in 



11] the other mechanism essentially contribute to this 



process practically all over the phase space. The thing is that in the leading order of pQCD 
the gluonic production described by 36 diagrams (see Fig. [3]) and only a small part of them 
belongs to the fragmentation type. The other diagrams can not be interpreted as b quark 
production followed by the fragmentation into B c meson (see for example the diagram (1) in 
Fig. [3]). One can say that these (recombination) diagrams correspond to the "independent" 
creation of cc and bb pairs. In spite of the suppression by a factor l/p% comparing with 
the fragmentation diagrams, their contribution into amplitude can not be neglected even at 
relatively large transverse momenta. As one can see in Fig. H] the fragmentation approach 
is valid only at transverse momenta larger than 5 -j- 6 masses of B c . This behavior does not 



contradict the factorization theorem 



151 ]. which predicts the fragmentation dominance at 



limit of infinite transverse momenta, but says nothing about transverse momentum value at 
which the nonfragmentation contribution becomes negligible. 



2 Within the valence quark approximation \B*) = {b^q^, (b^q^, + b_iq\)/^/2, b^q^} and \B) = ib^q^ — 
b_\.q\) / \/2. Thus, if it is supposed that all possible spin combinations are produced with the same proba- 
bility, then cr(B*)/a{B) = 3. 
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Figure 4: The cross section distribution on 
transverse momenta for the B c gluonic pro- 



Figure 5: The cross section value dependence 
on the interaction energy for the gluonic B c 



duction at interaction energy ^/s gg = 100 GeV production [8(]. 
obtained within pQCD 8] in comparison with 
the fragmentation model prediction (see for- 
mula ([5]) for the fragmentation model). 

The total gluonic cross section predicted using full set of leading order diagrams essentially 
differ from the fragmentation approach in absolute value as well as in shapes of interaction 
energy dependence (see Fig. [5]). It is worth to note that near the threshold the fragmentation 
prediction values are larger than the prediction values obtained within pQCD calculation 
because of the incorrect phase space counting in fragmentation approach. 

To obtain the cross section estimations for real hadronic experiments we need to convolute 
the partonic cross sections with the parton density functions: 

OhadronicO) = J J <Ty (s ij} Q R )fi{x h Q F )fj{Xj, Q F )dx l d Xj (8) 

As it is predicted within pQCD , about 90 % of the B c mesons at LHC energies 
will be produced in the gluonic fusion (Fig. |3]). Therefore in our pQCD calculations we can 
neglect the other partonic subprocess. 
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Figure 6: The cross section distribution over 
the B c transverse momentum for the process 
pp — > B c + X at y/s = 14 TeV within partonic 
approach ([8]) (see also 



Figure 7: The normalized cross section distri- 
butions over the invariant mass of B c (B*) and 
c-quark within fragmentation approach. 



The convolution of the gluonic subprocess cross section with the gluonic structure func- 
tions partially hides the differences between the pQCD predictions and the fragmentation 
approach. For example, one can see in Fig. |H]that at hadronic level the difference in shapes 
of px distributions is not so obvious, as at partonic level. 3 The predicted difference of the 
hadronic cross section values is about of 2. Obviously, such a difference is not essential for 
the calculations of forth order on a s . Nevertheless, the relative yield of B* and B c does not 
depend on a s and could indicate the production mechanism. Even in the kinematical region 
where the fragmentation model could be applied the value of -Rhadr = o"hadr(-B*)/ cr hadr(-Bc) 
predicted within pQCD is about 2.6 instead of 1.4 obtained within the fragmentation ap- 
proach. To measure this value one need to detect B* meson which with unit probability 
decays into B c and photon. However, it is quite difficult to detect such a process experi- 



3 To estimate the B c hadronic production cross section the CTEQ6L parametrization of gluonic structure 
function has been used on the scale Qf = 10 GeV 16 1. 



mentally due to the expected small mass difference between B* and B c mesons 
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AM = M B * - M Bc = 65 ± 15 MeV. (9) 
Therefore in laboratory system the maximum energy of emitted photon is 

Umax = (7 + V(l 2 ~ !)) 

where 7 is B* 7 factor. Even for B* with energy ~ 30 GeV u max ~ 0.7 GeV. Thus one can 
conclude that there is no certainty that the method based on separation of B* from B c can 
be used to study the B c production mechanism. This is why in the next chapter we suggest 
another way to research the B c production mechanism. 



III. THE CROSS SECTION DISTRIBUTION ON THE INVARIANT MASS OF 

D AND B c MESONS 

Within fragmentation approach the shape of the cross section distribution over the in- 
variant mass of B c and c-quark is roughly determined by 6-quark virtuality (see diagram 
(1) in Fig. [1]). This is why the distribution should be relatively narrow. Our analytical 
calculations confirm this supposition. Let us denote by di_+ Bc (M Bc+5 ) and di^ B *(M Bc+s ) 
the normalized cross section distributions over the invariant mass of B c (B*) and c-quark 
within fragmentation approach: 

d b ^ Bc (M Bc+5 ) = —-7^, (10) 



1 



da%^ 



d^B*(M M ) = (11) 



It can be easy to obtain from 12J that: 



a b i dM Bi 



mz. 



l B c +c 

c J V 

\2„ H/T2 



[1 - z)(l + rzfrM 2 B (2(1 - 2r) - (3 - 4r + Ar 2 )z + (1 - r)(l - 2r)z 2 )rM% c 



1 - (1 - r)z)\M Bc+ - c - m\f (1 - (1 - r)z)(Ml + - c - 

4r 2 (l - r)M Bc N 



(12) 
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c+c) 



8a s \R(0)\ 2 M Bt 

2771171% 



+ ' 'dze{M 2 Bc+5 



Ml 



l-z 



(1 -z)(l + 2rz + (2 + r 2 )z 2 )rM 2 Bc (2(1 + 2r) - (1 + Ylr - Ar 2 )z - (1 - r)(l + 2r)z 2 )rM 4 c 
(l-(l-r)^) 2 (M| c+c --m2) 2 " (l-(l-r)^)(M 2 * 



(Ml c+C - - m 2 ) 4 



l2r 2 (l-r)M% . 

1 , (13) 



where M is .B c meson mass, M Bc+5 is the invariant mass of B c (B*) and c-quark. 
Integrating over z leads to the following expressions: 

, v Aa s \R(Q)\ 2 M 2 Bc r(l + r) 2 M Bc+5 
d- b ^ Bc {M Bc+5 ) = ^ x 



277rm3(l -r) A A 2 



9 4Af£ r(l + r)(-l + r) 
1 + 2r 2 + ^ 



21 i f A + B , 



(1 + r) 2 



-{il4 c (-l + r) 8 (l + r) - Ml c (-1 + r) 4 (4 - 14r + 7r 2 + r 3 )Ml c+c - 



- M 2 c (-1 + r) 3 (l - 8r + r 2 )M 4 c+E + (2 + 4r - r 2 + r 3 )M Bc+£ } ■ B 
- 2M Bc+ - c {M Bc (-l + r) 4 (l - 4r + 2r 2 ) 

2M 2 c (-l + r) 2 (l - 2r + 2r 2 )Mi c+c - + (1 + 2r 2 )M 4 c+5 } • log {^T^ 



(14) 



. ,_ 4a s |i?(0)| 2 r(l + r) 2 M 2 e M Be+5 ^ 

dl^B t {M Bc+5 ) - 277rm 3 (1 _ r)M 2 X 



(-2^ t^K" 1 + r ) 2 ( 3 - 8r + 2r2 ) " (3 + 4r + 2r 2 )Ml c+c -] • log {^^} 



+ 



{M Bc (-l + r) 8 (l + r) - M Bc (-l + r) 4 (2 - 9r 2 + r 3 )M 2 Bi 



+c 



(1 + r) 2 

- M 2 c (-1 + r) 3 (-ll + 8r + r 2 )M 4 c+c - + (12 + 6r - r 2 + r 3 )M| c+c _} . B 
- 2M Bc+5 {M Bc (-l + r) 4 (3 - 8r + 2r 2 ) 



-2M 2 f- 



(-1 + r) 2 (3 -2r + 2r 2 )M 2 Bc+5 + (3 + 4r + 2r 2 )M 4 c+c -} • log j 



( A-B 



where 



and 



B 



A = {l+r){M 2 Bc+ - c -{l-r) 2 M 2 B ^ 
(1 " r) sj (M 2 c+5 - (1 - r) 2 M 2 c ) (m 2 c+c - - (1 + r) 2 M 2 c ) . 



(15) 

(16) 
(17) 
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Figure 8: The cross section distributions over the invariant mass of B c and D mesons calculated 
within pQCD for the process pp — > B c + X at = 14 TeV: without cuts (a) and for y/§ > 40 GeV 
(b). Solid curves: D c ^d = 5(1 — z); dashed curves: D C -^D = z 2 - 2 (l — z). 

The variation of di^ Bc and <%_> B , with Mb c +c are shown in Fig. [7J 

Here we face the problem how to transform the invariant mass of B c (B*) and c-quark 
Mb c +c to invariant mass of B c (B*) and D-meson M Bc+ f). Within fragmentation mechanism 
it naturally to assume that the D meson takes away the total momentum of c-quark, because 
the production of c quark is the last step of emission process. Such an assumption is not 
obvious for the recombination contribution. This is why two hadronization models of c quark 
have been chosen: 

1. D meson takes away the total momentum of c-quark: D c ^d = 5(1 — z); 

2. D meson takes away part of c-quark momentum according to Kartvelishviiy-Petrov- 
Likhoded fragmentation function: D c ^d = z 2 ' 2 (l — z) \l\ - 

The cross section distributions over M Bc+D . are shown in Fig. [S] for the different kinemat- 
ical regions. One can see that for the total phase space the cross section distribution looks 
like the obtained within the fragmentation approach. Nevertheless it become essentially 
wider for large interaction energies, whereas within the fragmentation approach it should 
remain the same. Thus one can conclude that the recombination contribution is dominant. 
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Figure 9: The cross section distribution over 
the cosine of angle between the directions of 
motion of B c and D mesons predicted within 
pQCD for the process pp — > B c + X at ^/s = 
14 TeV. 
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Figure 10: The dependencies of averaged c 
quark energy on B c meson energy are repre- 
sented for the different cuts on 6-jet transverse 
momenta. 



IV. THE ANGLE CORRELATIONS AND THE DECAY LENGTH RATIO IN 
THE ASSOCIATIVE PRODUCTION OF B c AND D MESONS. 

As it is shown in Fig. HI the cross section distribution on cosine of the angle between the 
B c meson and D meson has the sharp maximum at cos# ~ 1. Moreover, about a half of B c 
mesons is associated by the D meson moving in the close direction: 9 < 26°. Therefore, D 
meson can be used to detect B c meson. 

In the Fig. [10] the dependencies of averaged c quark energy on B c meson energy are 
represented for the different cuts on 6-jet transverse momenta. It can be concluded from 
these plots that at any 6-jet transverse momentum 

(Ec) > l-2E Bc . (18) 



D meson we obtain that 



(E D ) > 1.2E Bc 



(19) 
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for D c ^d = 5(1 — z) or 

(E D ) > 0.7E Bc (20) 

for D c _ D = z 2 - 2 (l-z). 

The decay lengths depend on particle energies and lifetimes as follows: 

(h) ^ -ct d l Bc ~ — -ct d . 21 

m D m Bc 

Taking into account that t d /t Bc ~ 2 we obtain: 

M > 5. (22) 
This value should be compared with the fragmentation model prediction: 



^frag 



1 -r- 2. (23) 



It worth to note that the ratio {Id)/Ib c can not be used to detect the B c meson because 
the distributions on Id at fixed l Bc very wide. Nevertheless the experimental measuring of 
this ratio could indicate the production mechanism. 



V. CONCLUSIONS 



The following conclusion can be drawn from the performed calculations: 

1. The cross section distribution over the invariant mass of B c and D meson depends 
essentially on kinematical cuts and can be used to research B c production mechanism 
at LHC. 

2. In many cases the B c and D mesons move in close directions. It could be useful to 
detect B c meson. 

3. The energies of B c and D mesons are comparable. The decay length of D meson 
by more than 5 times larger than the decay length of B c meson. The experimental 
research of the ratio between B c meson and D meson energies could shed light on B c 
production mechanism. 
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